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Abstract 

We conjecture that the highest energy cosmic rays, E > Eqzk, where Eqzk ~ 5 • 10^^ 
eV is the Greisen-Zatsepin-Kuzmin cut-off energy of cosmic ray spectrum, may provide a 
unique window into the very early epoch of the Universe, namely, that of reheating after 
inflation, provided these cosmic rays are due to decays of parent superheavy long-living 
X-particles. 

These particles may constitute a considerable fraction of cold dark matter in the 
Universe. We argue that the unconventionally long lifetime of the superheavy particles, 
which should be in the range of 10^'' — 10^'^ years, might require novel particle physics 
mechanisms of their decays, such as instantons. We propose a toy model illustrating the 
instanton scenario. 

Generic expected features of ultra-high energy extensive air showers in our scenario 
are similar to those of other top-down scenarios. However, some properties of the upper 
part of the cosmic ray spectrum make the instanton scenario distinguishable, at least in 
principle, from other ones. 



1. Among the top-down mechanisms attempting to explain the observations of uhra- 
high energy (UHE) cosmic rays beyond the Greisen-Zatsepin-Kuzmin cut-off , decays of 
primordial heavy particles are most obvious possibility. It is clearly more conventional — at 
least at first sight — than scenarios invoking topological or non-topological defects 0], though 
the latter may be a viable alternative. Heavy particles with lifetime of the order of the age of 
the Universe or greater may constitute (a substantial fraction of) cold dark matter, so the two 
observed features of the Universe — ultra-high energy cosmic rays and dark matter — may be 
related to each other. 

To get an idea of the range of properties of decaying particles (X-particles) that supposedly 
produce UHE cosmic rays, let us make the following simple observations. 

First, assuming sizeable hadronic component (jets) among the decay products, the flux of 
protons or gammas of energy E on the Earth is estimated as 

dlnE 47rrx ^'^ ^' d InE ' ^ ' 

where Nj is the number of jets in a typical decay, Rp^^ is the effective distance to X-particles, 
nx is the number density of X-particles at the scale Rp^^ , tx is the X-particle lifetime and 
dN / d InE is the fragmentation function. In the following estimates we take Nj ~ 1 — 10 (we will 
soon see that large jet multiplicity may be favored in some decay scenarios) and dNp^^/d InE ~ 
(a few) • (10 — 100) in the energy range of interest, E > few) • lO^'' GeV (the latter estimate 
comes from bold extrapolation of the fragmentation functions of ref.[0 to extremely high jet 
energies). For the effective distance we take R ^ 100 Mpc, with understanding that the 
actual value of R may be much smaller than 100 Mpc if X-particles are clumped. In fact, our 
conclusions will be fairly insensitive to the actual values of the above parameters. 
The second relation is 

mx <nx >= ^xpcrit, (2) 

where < nx > is the average number density of X-particles and Qx ^1- If X-particles 
are clumped, the density nx entering Eq.(^ may be several orders of magnitude larger than 
< ux >■ Again, this uncertainty will not affect our main conclusions, and we set ux ~ < nx > 
in what follows. In order to produce cosmic rays of energies E > {a few) • 10^^ GeV, the mass 
of X-particles is to be very large, mx ~ 10^'^ GeV. 

Let us now estimate the range of X-particle densities required. From Eq.(|^) we find a bound 
for the present X-particle density-to-entropy ratio 

nx/s < 10-^\ (3) 

On the other hand, to produce the observed flux of UHE cosmic rays, the density of X-particles 
should not be too small. Keeping in mind that their lifetime, tx, cannot be much smaller than 
the age of the Universe, 

r > 10^° yr, (4) 

we obtain from Eq.(^) 

nx/s > 10-33. (5) 

Even though the window for nx is very wide, the estimates (^ and (|^) raise the issue of the 
production of X-particles in the early Universe. 
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Alternatively, Eqs.(^ and may be used to place an upper bound on the lifetime of 
X-particles, 



Again, the window for tx is wide, but the estimates and indicate another problem, 
namely, that of the particle physics mechanism responsible for long but finite lifetime of very 
heavy particles. 

In the rest of this paper we propose possible scenarios for i) generating the abundance of 
X-particles in the range (1) by processes in the early Universe, and ii) explaining the 
lifetime of X-particles in the range (^, 

Before coming to our main points, let us stress that these two problems inherent in theories 
with very heavy and almost stable particles were realized long ago (see, e.g., ref.[^ and refer- 
ences therein) in different contexts. In ref. |^ it was proposed that the problem i) may be solved 
by large entropy generation in the Universe after the heavy particles freeze out of thermal equi- 
librium, while their long lifetime may be due to very large dimension of operators responsible 
for the decay. We leave for the reader to judge how exotic are alternative possibilities that we 
discuss below. 

2. If the temperature in the early Universe at some epoch exceeded the mass of X-particles 
and then decreased smoothly without large entropy generation, the freeze-out density of X- 
particles would greatly exceed the bound (0). A way out of this problem is provided by inflation 
and subsequent reheating. To explain the small abundance of X-particles, the reheating temper- 
ature Tr must be much smaller than mx-, so that X-particles were never at thermal equilibrium 
after inflation. In that case X-particles-to-entropy ratio is exponentially small. 



where the constant depends on a number of factors (the coupling constant responsible for pair 
production of X-particles, the effective number of degrees of freedom, the ratios mx/T^ and 
Mpi/iTLxi etc.) and is of order 10~^ with several orders of magnitude uncertainty. As the 
dominant suppression comes from the exponential factor, the reheating temperature can be 
estimated from (|), (|]) fairly precisely. 



and should be in the range 10^^ — 10^^ GeV, depending on mx- Note that this range is realistic 
in many scenarios of inflation. 

Hence, inflationary scenario can easily explain small value of the present density of X- 
particles in the space. Conversely, the determination of mx from measurements of the upper 
end of cosmic ray spectrum would allow for rather precise estimate of the reheating temperature. 
Ultra-high energy cosmic rays may indeed serve as a window to reheating epoch in the early 
Universe. 

3. Explaining long lifetime of X-particles is much harder. Conventional perturbative mech- 
anisms cannot be responsible for cosmologically large tx (unless very high dimension operators 
are involved 0), so one turns naturally to non-perturbative phenomena. A well known example 
is instantons that produce exponentially small effects in weakly coupled theories. If instantons 
are responsible for X-particle decays, the lifetime is roughly estimated as 



Tx ^ W yr. 



(6) 



nx/s = const-exp(— 2mx/Tr), 



(7) 




(8) 



Tx ^ mx ■ exp(47r/ax), 



(9) 
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where ax is the couphng constant of the relevant (spontaneously broken) gauge symmetry. 
From Eqs.(^), @ we find that the coupling constant (at the scale mx) is 



ax = . (10) 

^ 10 12 ^ ' 

Hence, we are lead to introduce additional non-Abelian gauge interactions with fairly large 
coupling constant at high energy scale. 

To illustrate this possibility let us consider a toy model with SU{^)x gauge interactions 
added to the Standard Model. The SU{2)x gauge symmetry is assumed to be broken at 
sufficiently high energy scale. Some conventional quarks and leptons carry non-trivial SU {2)x 
quantum numbers (say, SU (2)x may be right-handed subgroup of a left-right symmetric theory, 
or it may be a horizontal group with generations forming SU{2)x triplets). Let there be twoQ 
left-handed SU{2)x fermionic doublets X and Y and four right-handed singlets, all of which 
are singlets under SU{2)l x SU{?>)c of the Standard Model. After SU{2)x breaks down, all X- 
and y-particles acquire large masses in a manner similar to the Standard Model. We further 
assume that X and Y carry different global quantum numbers, so there is no mixing between 
them. 

Under these assumptions the lightest of X-particles and the lightest of F-particles (we call 
them X and Y at certain risk of confusing notations) are perturbatively stable. However, 
SU{2)x instantons induce effective interactions violating global quantum numbers of X and 
Y . Say, if X is heavier than F, then S'[/(2)x-instantons lead to the decay 

X ^Y + quarks + leptons (11) 

with the rate estimate given by Eq.(^. It is this type of processes that may be responsible for 
the production of ultra-high energy cosmic rays. 

Let us point out a few features of this scenario which seem generic. 

i) Decays induced by instantons typically lead to multiparticle final states. The number 
of quarks (jets) produced in the process (0) should be rather large, of order 10, and their 
distribution in energy should be fairly flat. In principle, the spectrum of cosmic rays within 
this scenario should be distinguishable from the spectrum predicted by other mechanisms (like 
two-jet decays of heavy particles born in the interactions of topological defects). Also, there 
are necessarily hard leptons among the decay products in the process (0). 

ii) If y-particles are indeed perturbatively stable, they are also stable against instanton- 
induced interactions (because of energy conservation and instanton selection rules). Then the 
dark matter in the Universe may consist predominantly of F-particles, while the admixture 
of X-particles is small. Since the abundance of F-particles is given by the formula similar to 
Eq.(^, and since the density of X-particles is bounded from above, Eq.(||), the mass splitting 
between X- and F-particles should not be large, 

nix < 2my. (12) 

Alternatively, the Higgs sector and its interactions with fermions may be organized in such a way 
that F-particles are in fact perturbatively unstable (while X-particles remain perturbatively 
stable). In that case the heavy candidates for dark matter are X-particles, and the approximate 
degeneracy (O) need not hold. 



^The SU{2)x anomaly prevents the number of SU{2)x doublets from being odd. 
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iii) Because of the instanton selection rules, this scenario for the slow decay of heavy par- 
ticles is rather restrictive. For example, X-particles cannot be colored (otherwise the SU (2)x 
instanton vertex would include at least three X-fields); if F-particles are stable, they cannot 
be colored either. This fits nicely to the expectation that dark matter particles do not expe- 
rience strong interactions. On the other hand, X-particles (and F-particles, if stable) cannot 
be weak doublets for the same reason, so the heavy dark matter in this scenario does not have 
electroweak interactions, too. 

To conclude, an explanation of ultra-high energy cosmic ray events beyond the GZK cut- 
off by decays of hypothetical heavy particles of cosmologically long lifetime is not unrealistic 
from both cosmological and particle physics points of view. Detailed study of the upper end 
of cosmic ray spectrum will provide insight into the decay mechanism involved, and allow for 
the determination of the mass of X-particles. As the latter has been argued to be related 
to the reheating temperature, the ultra-high energy cosmic rays may become a clue to the 
end-of-inflation epoch in the early Universe. 

After this work was presented at this Conference, we received a paper by Berezinsky, Kachel- 
riess and Vilenkin ^ where the decays of heavy particles have been also considered as the origin 
of cosmic rays beyond the GZK cut-off. Their main point is that X-particles are expected to 
concentrate in the galactic halo, and in this way one easily avoids the constraints coming from 
the analysis of the cascade radiation 0. Their proposal for long but finite lifetime of X-particles 
is that it is due to quantum gravity (wormhole) effects. 

We are indebted to V. Berezinsky, G. Farrar, W. Ochs, G. Pivovarov, S. Sarkar, D. Semikoz, 
G. Sigl and L. Stodolsky for stimulating discussions. The work of V.K. was supported in part by 
Russian Foundation for Basic Research grant 95-02-04911a. The work of V.R. was supported in 
part by Russian Foundation for Basic Research grant 96-02-17449a and U.S. Civilian Research 
and Development Foundation for Independent States of FSU (CRDF) award RPl-187. 
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